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\ We present a theoretical formalism by which the global and the local mass 

j functions of dark matter substructures (dark subhalos) can be analytically es- 

timated. The global subhalo mass function is defined to give the total number 
density of dark subhalos in the universe as a function of mass, while the lo- 
C*~) \ col subhalo mass function counts only those subhalos included in one individual 

\ host halo. We develop our formalism by modifying the Press-Schechter theory 

to incorporate the followings: (i) the internal structure of dark halos; (ii) the 
{T} ', correlations between the halos and the subhalos; (iii) the subhalo mass-loss ef- 

fect driven by the tidal forces. We find that the resulting (cumulative) subhalo 
mass function is close to a power law with the slope of ~ — 1, that the subhalos 
contribute approximately 10% of the total mass, and that the tidal stripping 
effect changes the subhalo mass function self-similarly, all consistent with recent 
numerical detections. 



Subject headings: cosmology: theory — large-scale structure of universe 



1. INTRODUCTION 



The dark halo substructures (dark subhalos) are the dynamically distinct, self-bound 
objects in virialized dark matter halos. The presence of substructures in the dark matter 
halos is a generic picture of the cold dark matter (CDM) cosmology. Recent numerical 
simulations of ultra-high resolution indeed confirmed that the dark halos are not smooth 
structureless objects but clumpy systems marked by a wealth of substructures (Tormen et 
al. 1998; Klypin et al. 1999; Okamoto & Habe 1999; Ghigna et al. 2000; Springel et al. 2002; 
Zhang et al. 2002; De Lucia et al. 2003; Hayashi et al. 2003; Zenter & Bullock 2003). 

Recently, the mass function of dark subhalos has drawn sharp attentions (Fujita et al. 
2002; Sheth 2003; Blanton 2003) especially because of its connection to the galaxy luminosity 
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function. Yet, it is not an easy task to derive the subhalo mass function either in numerical 
or analytical ways. The numerical approach to the subhalo mass function using N-body 
simulations still suffers from resolution effects related to the so called over-merging problem 
(Klypin et al. 1999). Even recently available ultra-high resolution simulations are capable 
of producing only the local subhalo mass function, i.e., the mass function of the subhalos 
within one individual dark halo (Okamoto & Habe 1999; Ghigna et al. 2000). Given the 
importance of the subhalo mass function as a clue to understanding of the galaxy luminosity 
function, however, what is also desired is the global subhalo mass function, i.e., the mass 
function of all the subhalos in the universe, irrespective of the host halos. 

As for the analytic approach, the hindrance is the complexity of the subhalo evolution. 
For the mass function of dark halos, we already have a remarkably successful theory devel- 
oped by Press & Schechter (1974, hereafter PS). The principle of the PS theory is this: the 
formation and evolution of dark halos can be traced by the linear theory, assuming (i) dark 
halos have no internal structure; (ii) dark halos form independently of their surroundings; 
(hi) dark halos do not lose mass in the evolution but only hierarchically merge via gravity. 
Unlike the case of the halo mass function, however, the subhalo mass function cannot be 
derived under such simple assumptions. The subhalos are, by definitions, the internal struc- 
tures of the halos, being placed in highly dense surroundings, and thus the formation and 
evolution of the subhalos must depend strongly on their surroundings. Among the various 
consequences from the surrounding influences, the most significant one is the subhalo mass- 
loss: the subhalos do not only gravitationally merge but also get disrupted or at least lose 
considerable amount of their mass through the interaction with the surroundings. In fact, it 
has been demonstrated by several N-body simulations that the subhalos lose most of their 
mass throughout the evolution, contributing after all only 10 — 15% of the total mass of the 
host halos (e.g., Tormen et al. 1998). There are three different processes that can drive the 
subhalos to lose mass: the global tides generated by the host halos, the dynamical frictions, 
and the close-encounters with the other subhalos. Apparently, the subhalo mass-loss is quite 
a complicated process, so that it would be practically impossible to take into account its 
effect fully in deriving the subhalo mass function analytically. That was why all previous 
analytic approaches had to make the unrealistic assumption that the subhalos do not lose 
mass during the evolution (Fujita et al. 2002; Sheth 2003; Blanton 2003). 

However, the mass-loss phenomenon is the most essential feature of the subhalo evo- 
lution, which must be taken into account in order to estimate the subhalo mass function 
in any realistic sense. Here we attempt for the first time to estimate both the global and 
the local subhalo mass functions with the subhalo mass-loss effect taken into account. To 
make the theory analytically tractable, we still make some simplified assumptions that the 
subhalo mass-loss is mainly driven by the global tides, and that the condition for a subhalo 
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to survive the global tides is a simple function of the distance from its host halo. 



2. 



FORMALISM 



The (differential) global subhalo mass function, dll | M d\nM s , is defined as the number 
density of subhalos in logarithmic mass range [lnM s ,lnM s + d\nM s ]. To estimate it, we 
assume the following. 

(i) The gravitational collapse process to form dark matter halos follows the Top-hat spher- 
ical dynamics, according to which a dark halo of mass M forms at redshift z if the 
density contrast 5 (5 = Ap/p, p: the mean mass density of the universe) of a La- 
grangian region in the linear density field smoothed on mass scale of M satisfies the 
gravitational collapse condition of 5 = S c (z) where S c (z) is the critical density contrast 
at redshift z whose value depends on cosmology (Kitayama & Suto 1996). For a flat 
universe of closure density, 5 c (z) « 1.68(1 + z) (Gunn & Gott 1972). 

(ii) The mass function of dark halos is well evaluated by the PS theory as 



where v = 5 c /a(M) and a(M) is the rms density fluctuation of the linear density field 
on mass scale of M. 

(iii) A dark halo hosts a multiple of subhalos, each of which rotates on a circular orbit. 
The global tidal field of the host halo strips the subhalos, which drives the subhalos to 
either get completely destroyed or survive but with reduced mass. A subhalo of initial 
mass Mi before the tidal stripping effect at an initial distance r from its host halo of 
mass M 2 eventually survives the tidal stripping effect, ending up with having reduced 
mass of M s (< Mi), if Mi, M 2 , and r, satisfy the condition 



where the proportionality constant c t is a free parameter. 

(iv) The spatial distribution of the subhalos inside their host halo follows that of the dark 
matter particles, i.e., the profile given by Navarro, Frenk, & White (1996) (hereafter, 
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NFW): 



P M2 (r) oc 
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(r/r s )(l + r/r s y 



where r s is the scale radius. 
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It is worth noting that equation (2) is reminiscent of the familiar tidal-limit approxima- 
tion, according to which a subhalo in the tidal field loses all mass beyond its tidal radius, r t . 
We set c t as a free parameter since its precise value depends on underlying assumptions: if 
the halo potentials can be approximated as point mass and the linear size of the subhalo is 
much smaller than the distance from the host halo, the tidal radius has a simple expression 
of r t = [ctMxjM^l^r with Ct — 1/2 (the Roche limit); if the effect of the centrifugal force is 
taken into account, then c t = 1/3 (the Jacobi limit); if the halos are treated more realistically 
as the extended mass profiles, then the tidal radius has a more general expression (see, e.g., 
Tormen et al. 1998). Although equation (2) is an obvious oversimplification of real tidal 
mass-loss process (Hayashi et al. 2003), there are numerical clues that the fraction of the 
survival subhalos has a strong correlation with the distance from the host halos (Okamoto 
& Habe 1999). It implies that the tidal survival condition should depend on the subhalo 
orbital distance as well as the host halo and the subhalo mass. Therefore, equation (2) may 
be the simplest possible choice for the functional form of the tidal survival condition. 

Let us first consider Pm 2 (M\, r), the conditional probability that a subhalo has an initial 
mass greater than Mi (before the tidal mass-loss) provided that it rotates upon a host halo 
of mass M 2 at a distance r. According to the hypothesis (i), Pm 2 (^i! r ) — Ps 2 =s c [Si(r) > S c ] 
where 5i (with i — 1,2) is the density contrast of a Lagrangian region in the linear density 
field smoothed on mass scale M i: and S ci is the critical value of Si. The conditional probability 
Ps 2 =s c [Si(r) > S c ] can be computed from the Gaussian probability density distribution with 
the help of the Bayes theorem. For the case of the sharp k-space filter, it has the following 
simple analytic form (Yano et al. 1996): 



1 f°° f x 2 \ 
P M2 (M 1 ;r) = —j= J exp ^ — — J dx, with (3 = 
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Here of (with i = 1,2) and cr 2 (r) are the mass variance of the linear density field on mass 
scale of Mi, and the linear density cross correlation, respectively, given as 

A\k)d\nk, a 2 c = A\k)^^d\nk (5) 

where is the dimensionless power spectrum of the linear density field, 7 = a 2 / (17^2), 
and the integral upper limit k C i is related to Mj by k C i = (67r 2 p/Mj) 1//3 . In fact, it was 
Yano et al. (1996) who first incorporated the spatial correlations between the dark halos 
themselves into the PS theory. Fujita et al. (2002) used the formalism of Yano et al. (1996) 
to estimate the local subhalo mass function with the spatial correlations between the host 
halos and the subhalos taken into account. However, Fujita et al. (2002) assumed that the 
spatialdistribution of the subhalos is uniform, and averaged Pm 2 (Mi, r) over r without taking 
into account the tidal mass-loss and its correlation with r. 
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Next, let us consider P M2 {r), the probability of finding a subhalo in a spherical shell of 
radius r and thickness dr around a host halo of mass M 2 in the Lagrangian space. According 
to the hypothesis (iv), it can be written as 



M 2 



(r/r s )(l + r/r s ) 



7 4:7ir 2 dr, with A 



47rr3[ln(l + £ ) 



R2 



(6) 



r s +R 2 J 



where i? 2 is the virial radius of M 2 , and the amplitude A is determined to satisfy the nor- 
malization constraint of j^ 2 PM 2 ( r ) — 1- We approximate _R 2 by the Top-hat radius of M 2 , 
and adopt the empirical relation for r s proposed by Klypin et al. (1999): 



3M 2 
4irp 



nl/3 



124 



h- x M t 







0.084 



(7) 



It is also worth mentioning here that 5,R 2 , r, and r s are all measured in the Lagrangian 
space where the density field is still Gaussian. 

The joint conditional distribution, Pm 2 {j-, Mi), i.e., the probability of finding a sub- 
halo of mass greater than M\ at a distance r from a host halo of mass M 2 , can be derived 
from equations (4) and (6) by using the Bayes theorem: Pm 2 {Mi,t) = Pm 2 {Mi, r )PM 2 { r )- 
The partial derivative, dPM 2 {r, M 1 )/9M 1 , is proportional to the fraction of the host halo 
volume occupied by those subhalos of mass Mi at a distance r, f'M 2 {Mi,r), such that 
fM 2 (Mi,r) = (M 2 /p)\dPM 2 (r,M 1 )/dMi\ where the proportionality factor (M 2 /p) is noth- 
ing but the average volume of the host halo. Following the familiar PS-like approach, the 
initial number density of the subhalos of mass Mi at a distance r inside a host halo of 
mass M 2 equals f M2 (Mi,r) divided by the average volume of the subhalo Mi/p such that 
dNM 2 {r, Mi)/(47rr 2 dr dM{) = 2 (p/Mi) fM 2 {Mi) where the factor 2 is the normalization con- 
stant introduced by PS (see also, Peacock & Heavens 1990; Bond et al. 1991; Jedamzik 1995; 
Yano et al. 1996; Fujita et al. 2002). We end up with the following expression: 



d 2 N M2 (r,Mi) 
Anr 2 dr dlnMi 



1 

Wi 



rflnai 



din Mi 



1 - 7 2 / (r/r s )(l + r/r s )' 



exp - 



e 



(8) 



where Mi = M x /M 2 and f = 



1 — 7^7) with Vi = bdj a i for % — 1,2. Equation (8) is 



the conditional mass and spatial distribution of the subhalos provided that they are included 
in the host halos of mass M 2 before the tidal mass-loss. If the subhalo mass were conserved, 
the local subhalo mass function would be computed simply by integrating equation (8) over 
r. However, the subhalos in reality either get destroyed or lose mass, so that the subhalo 
mass function cannot be simply obtained in that way. Notwithstanding, the reduced mass 
of survived subhalos can be determined from the informations of Mi, M 2 , and r according 
to the hypothesis (iii). 
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Hence, we find the local subhalo mass function after the tidal mass loss as 




where 5p represents the Dirac-delta function, and r c represents the lower limit for the subhalo 
survival: if a subhalo is located at a distance smaller than r c , they get completely destroyed 
by the strong tidal stripping effect. The value of r c has been empirically found to be a few 
times the scale radius r s (Ghigna et al. 2000; Hayashi et al. 2003). The upper limit R2 
in the integration of r is set from the expectation that the subhalos should be inside the 
virial radius of the host halo. Note that equations (8) and (9) are all expressed in terms of 
the rescaled subhalo mass Mi = Mi/M 2 . Their dependence on the host halo mass comes 
only implicitly from 7 and £. It implies that the subhalo mass function should not depend 
strongly on the host halo mass, consistent with numerical detections (e.g., De Lucia et al. 
2003). Furthermore, we compute the contribution of the subhalos to the host halo mass by 
integrating equation (9) over M s , and find that it is approximately 10% of the host halo 
mass, which is also consistent with numerical simulations (e.g., Tormen et al. 1998). Finally, 
the global subhalo mass function is obtained by multiplying the local subhalo mass function 
by the host halo mass function and integrating it over the host halo mass: 



where the host halo mass function dN "(M2) / 'din M2 is given in equation (1). 

Figure 1 plots the cumulative local (upper panel) and the global (lower panel) mass 
functions at redshift z = for the case of for a flat CDM cosmology (Bond & Efstathiou 
1984) with the following choice of the cosmological parameters: the matter density fl m = 0.3, 
the vacuum energy density Q\ = 0.7, the shape parameter T = 0.5, the dimensionless Hubble 
constant h = 0.5, and the rms density fluctuation on the scale of 8/i _1 Mpc, a s = 0.7. For 
the local subhalo mass function, the host halo mass is chosen to be M 2 = 10 14 M & h~ l . In 
each panel, the solid and the dashed lines correspond to the cases of q = 1/3 and c t = 1/2, 
respectively. For comparison, the subhalo mass function with no tidal effect is also plotted as 
dotted lines. It is clear from Figure 1 that both the global and the local cumulative subhalo 
mass functions are close to a power law N(M S ) ~ M l s , and that the tidal mass-loss effect 
changes the subhalo mass function in a self-similar manner. We find that the power-law slope 
of the local subhalo mass function is I ~ —0.8 while that of the global subhalo mass function 
is slightly sharper I ~ —0.9. The power-law slope of the global subhalo mass function gets 
sharper in the high mass section. It is due to that massive subhalos experience the stronger 
tidal stripping effect (eq. [2]) to lose more mass, and also that the dark halos which can 
afford to hosting such massive subhalos are rare whose number density decrease sharply (eq. 




(10) 



[I])- 
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3. DISCUSSIONS AND CONCLUSIONS 

We provided for the first time a theoretical formalism in which one can estimate an- 
alytically the global and the local mass distribution of dark matter subhalos that undergo 
tidal mass-loss process. Adopting the simple tidal-limit approximation, we showed that the 
resulting mass functions are consistent with what has been found in recent high-resolution 
N-body simulations, providing theoretical clues to the unique properties of the subhalo mass 
distribution: the power-law shape, weak dependence on host halo mass, self-similar change, 
and roughly 10% contribution of subhalos to the total mass. 

Yet, it is worth mentioning that our subhalo mass functions are subject to several 
caveats. The most obvious one is that we have oversimplified the subhao mass-loss process, 
using the simple tidal limit approximation, and also ignored the effects of dynamical frictions 
and close encounters between subhalos. It has been shown by numerical simulations that 
the tidal limit approximation underestimates the subhalo mass-loss considerably (Hayashi 
et al. 2003). Although it was shown by simulations that the most dominant force that leads 
to the mass loss of the subhalos is the global tides (Okamoto & Habe 1999), the dynamical 
frictions and subhalo close-encounters may change the subhalo orbits making the subhalos 
more susceptible to the tidal forces (Tormen et al. 1998). We have also assumed simply that 
the subhalos rotate on stable circular orbits. However, in reality, the the subhalo orbits are 
quite eccentric, changing with time (Tormen et al. 1998; Hayashi et al. 2003). Definitely, it 
will be quite necessary to refine our formalism by making more realistic treatments of the 
subhalo evolution, especially its mass-loss process. 

Finally, we conclude that our formalism is expected to provide an important first step 
toward realistic modeling of the abundance distribution of dark halo substructures. 

We are very grateful to the anonymous referee who helped us greatly improve the original 
manuscript. We also acknowledge gratefully the research grant of the JSPS (Japan Society of 
Promotion of Science) fellowship. This research was supported in part by the Grant-in-Aid 
for Scientific Research of JSPS (12640231). 
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Fig. I. — The mass functions of the dark halo substructures at redshift z=0: Upper Panel: 
the local distribution and Lower Panel: the global distribution. The solid and dashed lines 
correspond to the two different values of the free parameter c t in the the tidal survival 
condition (see, eq. [2]). The dotted lines correspond to the case of no tidal stripping effect. 



